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ABSTRACT 

Recent timing observations of the double pulsar J0737— 3039A/B have shown that 
its transverse velocity is extremely low, only lOkm/s, and nearly in the Plane of 
the Galaxy. With this new information, we rigorously re-examine the history and 
formation of this system, determining estimates of the pre-supernova companion mass, 
supernova kick and misalignment angle between the pre- and post-supernova orbital 
planes. We find that the progenitor to the recently formed 'B' pulsar was probably 
less than 2 Mq, lending credence to suggestions that this object may not have formed 
in a normal supernova involving the collapse of an iron core. At the same time, the 
supernova kick was likely non-zero. A comparison to the history of the double-neutron- 
star binary B 1534-1- 12 suggests a range of possible parameters for the progenitors 
of these systems, which should be taken into account in future binary population 
syntheses and in predictions of the rate and spatial distribution of short gamma-ray 
burst events. 



1 INTRODUCTION 



. The double pulsar P SR J0737-3039A/B jBurgav et alJ 
■ |2003l:lLvne' et al .120041) is an outstanding laboratory for tests 
, of general relativity and of binary evolution and supernova 
1 theories. With its short orbital period, Pb = 2.4 hours, 
and moderate eccentricity, e — 0.088, the general-relativistic 
modifications to the Keplerian orbit are the largest known. 
Pulsars 'A' and 'B' have spin periods of 22.7 msec and 2.7 s, 
respectively; since both are observed as radio pulsars, the 
meiss ratio is obtained directly from pulsar timing, leading 
to new theory-independent constraints on strong-field grav- 
ity llLvne et al.ll2004l: iKramer et aDbOOeTl . 

Such double-neutron-star (DNS) systems are descended 
from high-mass X-ray binaries (HMXBs) in which both stars 
are massive enough to undergo supernova (SN) explos ions 
(e.g.. lTauris fc van den Heuvelll2006l; iDewi fc Po'3l2003l) . In 
brief, the J0737— 3039A/B binary is thought to have begun 
as two main-sequence stars with masses of at least 8 Mq . 
After a first mass transfer stage, the primary formed a neu- 
tron star in a core-collapse SN. The secondary evolved, and 
matter was accreted by the neutron star in an HMXB phase. 
Eventually, the secondary's envelope enlarged to meet the 
neutron star, which spiraled in, ejecting the secondary's en- 
velope. Angular momentum transferred to the neutron star 
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spun it up to a period of a few tens of milliseconds; it is now 
observed as the A pulsar. Its envelope expelled, the helium 
core of the secondary remained in a circular orbit around 
the neutron star until a second supernova left the B pulsar. 
The radio lifetime of the recycled A pulsar is far longer than 
that of the high-magnetic-field B pulsar; this is the reason 
most double-neutron star systems are observed with only 
the recycled pulsar still active. We are fortunate to observe 
the double pulsar during its relative youth. 



The observed binary orbital elements and the space mo- 
tion of the binary preserve important details of the evolu- 
tion, including the mass lost in the second supernova (and 
hence the mass of the exploding helium star) , the size of the 
pre-supernova orbit, and any asymmetry in the explosion 
itself. The lower limit on the mass of an exploding helium 
core is of particular intere st; evolutionary models sug gest it 
should be around 2.1 Mq llNomot Jl984l;lHabetJl986f) . Also 
of interest is the misalignment angle S between the first-born 
NS's spin axis and the post-SN orbital angular momentum, 
equivalent t o the tilt betwee n the pre- and post-SN orbital 
planes (e.g.- IWex et~ai]l2000h . Constraints on all of the evo- 
lutionary parameters are important for estimating the de- 
tectability of DNS systems in pulsar surveys, and hence for 
estimates of the total DNS birthrate and population. These 
quantities are of wide interest for stellar evolution and nucle- 
osynthesis, but particularly because of their importance to 
event rate estimates for gravitational- wave detectors (e.g.. 
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|Phinnevll99ll : iKaloeera et aljEoO^) and as possibl e progen- 
itors of short gamma-ray bursts fGRBs: IPacvnskil llQSGl. 

Early attempts to constrain the pre-SN parameters and 
the kick velocity of J0737— 3039A/B have led to ambigu- 
ous results. Shortly after th e discovery, analyses based on 
the orbital elements alon e IIDewi fc van den Heuvell 12004 
IWillems fc Kalogerdl2004D suggested that moderate kick ve- 
locities (of 60km/s or more) were needed, and that the 
pre-SN helium star was likely low-mass and overflowing its 
Roche Lobe at the time of the explosion. Later analyses by 
Willems et al. (2004, 2006) used observational constraints 
on the transverse velocity of the system, first ~ 140 km/s 
from scintillation measurements (iRansom et a l. 2004) and 
later < 30km/s from pulsar timing jKrameT et aL .20051) . 
to trace the motion of the binar y back through the gravi - 
tational potential of the Galaxy iKuiiken fc Gilmorelll989f) 
and derive probability density functions for the most likely 
kick velocit y, pre-SN mass, and tilt angle. In the most re- 
cent paper iWillems et al .l200d) . population synthesis mod- 
els were also used to estimate the most likely radial veloc- 
ity of the binary (which is not directly measurable); this 
helped to further restrict the allowed parameters. Overall, 
these authors favor the standard formation scenario, rea- 
sonably high kick velocities (70-180 km/s) and pre-SN mass 
m2.i of at least 2AIq, but acknowledge that for very small 
transverse velocities (~ 10 km/s) lower-mass progenitors are 
allowed or even favored. Piran fc Shaviv (2005a,b) have ar- 
gued, based on the assumption that the system will oscillate 
vertically about the Plane of the Galaxy, that its current lo- 
cation close to the Plane implies that it must have a low 
transverse velocity and have experienced a very small or no 
kick at birth, indicating a pre-formation B mass less than 
2 Mq and possibly an entirely different formation mecha- 
nism, s uch as the electron -capture supernova proposed by 
[Podsiadlowski et al.l (|2005h . 



2 NEW RESULTS 

We are now in a position to reconcile these conflicting re- 
sults. Our ongoing timing observations at the Parkes, Jodrell 
Bank Lovell and Green Bank Telescopes have yielded a well- 
measured proper motion, fj, — 4.2 ± 0.4 mas/year, directed 
toward a celestial position angle (North through East) of 
308.4 ± 6.5°, or a G alactic position angle of 247.9 ± 6.5° 
iKramer et alj2006h - that is, its transverse motion is nearly 
parallel to the Plane of the Galaxy. Our tentative detection 
of the timing parallax is consistent with the adopted dis- 
tance estimate_of_520pcd^ from the dispersion mea- 
sure llCordes fc Lazidl2002h . Our measurements thus yield 
an extremely low transverse velocity, only 10 km/s, allowing 
us to derive a better estimate of the history of the system's 
motion in the Galaxy. 

To determine not just the limits of the allowed progeni- 
tor mass and kick velocity, but also their most likely values, 
we adapt the an alysis that we previou sly used for the DNS 
PSR B1534-hl2 llThorsett et al.ll2005l. hereafter TDS05) to 
calculate posterior probability density functions (pdfs) for 
several choices of priors on observable or potentially observ- 
able parameters. 

The unknown physical parameters are: the pre- 
supernova companion mass m2,i, the pre-supernova or- 



bital separation at, the 3-dimensional supernova kick vector 
Vk, the 3-dimensional centre-of-mass velocity of the post- 
supernova system Vcm, the 'tilt' angle 5 between the pre- 
and post-supernova orbital planes, the angle fl of the cur- 
rent binary's line of nodes on the sky, and the sign of the co- 
sine of the current system's orbital inclination angle i (note 
that 0° < i < 180° and that sini is well measured through 
timing). Two components of the current 3-d velocity are 
now measured through timing. Two relations between m2,i, 
ai, and Vk can be obtai ned through conservation of en- 
ergy and mornentum (e.g. iKaloeeral W9t : IWex et al]|2000l : 

1 Willems et ai]|2006^ . lwillems et alJ (l2006^ derive the set of 
possible solutions for Vk, weighting by a likelihood that as- 
sumes a uniform prior on the magnitude 14 and an isotropic 
angular distribution. They then convolve this function with 
their pdfs for the system's radial velocity Vr, and proba- 
bility that a given system will move to the pulsar's current 
location. 

We formulate the problem slightly differently, using our 
insight from TDS05 that if Vr (and hence Vcm), 5, i} and 
cos i are measured or specifled, then at every potential birth- 
site in the Galaxy, the equations connecting the remaining 
variables reduce to a simple quadratic in m2,i, with at most 

2 possible real solutions. This allows us to assume prior dis- 
tributions for Vr and the potentially observable 5 and SI. 
Our 3-d pdf may be written: 

piVr,n,s\D,i) 

= p{Vr, SI, S, Vpcc, cos j|mi, 7712,/, a/, e/, Vt,I) 
oc p{mi,m2.f,af,ef, Vt\Vr, fl, S, Vpec, cos i, I) 
xp{Vr, r^, S, Vpcc, cos i| J) 

where mi is the mass of pulsar A, m2,/ is the mass of pulsar 
B, a/ is the post-SN semi-major axis, and e/ the post-SN 
orbital eccentricity. Vpec represents the 3-dimensional pe- 
culiar velocity of the system before the second explosion, 
and is treated as a nuisance parameter along with cosi. The 
likelihood p(mi, m2./, af,ef,Vt\Vr,^,5, Vpcc, cos i, I) is sim- 
ply 1 for each acceptable real positive-mass solution of the 
quadratic equation for that choice of parameters, and oth- 
erwise. 

Our analysis proceeds as follows: For each of a large 
number of trials, we pick a radial velocity Vr from a prior 
distribution. As discussed by Piran fc ShavivlteOO^ and rec- 
ognized by Willems et al. (2006), the choice of this prior has 
important effects on the final pdfs. For the purposes of illus- 
tration, we investigate two prior dist ributions: 1) gau ssian 
in Vr with a dispersion 200 km/s (Wi llems et alj|2006ll . and 
2) Vr = Vt/ ta.nO, where cos 6 is chosen uniformly between 
— 1 and -1-1. The second prior is the one preferred by Piran 
fc Shaviv and is arguably the more logical prior to use^. We 
consider the first prior to be extremely conservative, in that 
its dispersion is much larger than the observed transverse ve- 
locities of DNS systems. We sample the proper motion from 



^ The Willems et al. (2006) definition of C differs from ours, but 
as both studies assume uniform priors this is not very important. 
^ We note that the expected radial velocity of the pulsar's LSR 
to that of the Sun is very small, on the order of 5km/s. Given 
the uncertainties in, for example, the gravitational potential of 
the Galaxy, we consider it appropriate to discuss radial velocities 
centred on Okm/s. 
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gaussian distributions using the measured uncertainties on 
magnitude and direction, and the pulsar distance assuming 
a 20% gaussian uncertainty. 

We then follow the motion of the binary system back 
in time through the Galaxy, similarly to the procedure in 
IWillems et"airi|2n0fih and TDS05, incorporating the Sun's 
position and peculiar motion. We accept any position with 
\z\ < 50 pc and galactocentric radius R < 15kpc as a pos- 
sible birth site for the B pulsar. We integrate back in time 
the orbital ecce ntricity and s emi-major axis, according to 
the equations of IPetera I^G^), testing ages up to 100 Myr. 

At each birth site, 1000 trial sets of parameters 
are selected. The peculiar velocity Vpcc is assumed to 
be a Maxwellian with a 1-dimensional dispersion of 
12km/s, based on the proper motions of Be /X-ray binarie s 
JChevalier fc Ilovaiskvl Il998l: Ivan den Heuvel et al.1 l200Cl) . 
while the cos i < and cos i > cases are treated as equally 
likely. We draw each of fl and S from uniform distributions: 
0° < n < 360° and 0° < 5 < 180°. The angle 6 is poten- 
tially observable but currently unknown, although the lack 
of profile shape changes in A ^Manchest er et al1l2005ll im- 
plies that it is probably small. T he angle can in principle 
be estimated from scintillation dColes et alj 1200511 but its 
value depends strongly on the modeling of the anisotropy of 
the interstellar medium and is not currently well constrained 
(Coles & Rickett, private communication). For each set of 
parameters, we construct the quadratic equation and deter- 
mine whether there are solutions. We record progenitor and 
kick parameters and construct pdfs via histograms. 

An important point neglected by both Willems et al. 
and Piran & Shaviv is that a neutron star's mass con- 
tains a large negative contribution from gravitational self- 
energy, meaning that the minimum progenitor mass must 
be equal to the current neutron-star mass plus the magni- 
tude of the binding energy Eb- This quantity can be esti- 
mated, to within about 2 0%, as Eb ~ 0.084(M/Mp ) ^ Mq 
jLattimer fc YahiJ Il989l: iLattimer fc PrakashI 120011) . The 
majority of plausible equations of state predict a very 
slightly lower binding energy for a 1.25 Mq neutron star. 
The net implication is that a minimum (baryonic) mass of 
about 1.37 is required for the progenitor of B, and we 
consider only solutions which yield progenitor masses above 
this limit. 

To place the double pulsar in context, we compare it 
to the only other DNS with strong constraints on the pro- 
genitor, PSR B1534-fl2 (TDS05). (Recent proper motion 
measurements for PSR B1913-I-16 will n ecessitate a revis ed 
set of kick constraints for that pulsar (iMceeraOliQol.) 
For a fair comparison, we have revisited our analysis of 
B1534-M2, calculating pdfs as for J0737-3039A/B, while 
allowing for progenitor masses down to 1.48 Mp, and not us - 
ing the scintillation constraint on Q iBoedanov et al.ll2002ll 
in case anisotropy systematics also affect this measurement. 
Because of the strong misalignment angle constraint, fewer 
sets of trial parameters are tested per birth site, while more 
separate radial velocity trials are used. 



3 DISCUSSION 

Because of the constraint equations used, each point in the 
(Vk, rn2,i, ai) parameter space corresponds to a point in 
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Figure 1. Normalized posterior pdfs for DNS progenitors. In each 
panel, the red and black lines represent PSR J0737— 3039A/B 
progenitors, with the red line corresponding to the gaussian Vr 
prior with dispersion 200 km/s, and black to the Vr prior which 
is uniform in the cosine of the angle between the velocity vec- 
tor and our line of sight. The blue lines show pdfs for progeni- 
tors of PSR B1534+12, for the gaussian Vr prior with dispersion 
200 km/s. Only parameter ranges with significantly non-zero pdf 
values are shown. Note the association of large 5 with large kicks, 
as expected. 



the {Vr, 5, Q,) parameter space. Thus we have effectively 
constructed the posterior pdfs for the (Vk, ra2,i, ai) pa- 
rameters as well, and may marginalize to derive confidence 
ranges on all the parameters. For both J0737— 3039A/B 
and B1534-f 12, our pdfs for the physically interesting quan- 
tities are plotted in Figure Q and confidence ranges are 
summarized in Tabl e [J We recognize, as pointed out by 
I Willems et al.l ll2006l) . that projection effects are neglected 
in these 1-d pdfs. The most important such effect for 
J0737— 3039 is the association of large-Vk, large-m2,i and 
large-(5 solutions exclusively with the large Vr (^ 100 km/s) 
progenitors accessible only with the gaussian Vr prior. Be- 
cause these imply that the current velocity makes a very 
small angle to our line of sight, we regard them as inher- 
ently less plausible. 

We can draw several conclusions about the 
J0737— 3039B progenitor and the supernova kick. First, the 
explosion was asymmetric, but not extremely so. The Vk pdf 
extends down toward Okm/s for both Vr priors, but a non- 
zero kick is strongly favoured. That the supernova explosion 
was most likely asymmetric is bolstered by the B pul- 
sar's spin-orbit misalign ment observed throu gh long-term 
changes in the B profile jBureav et al."2005') and derived 
through modeling of the A eclipses (Lvutikov & Th ompsoiil 
l2005llLvutikwll2005ll . For the gaussian Vr prior, the kick 
tends to be directed out of the plane of the pre-SN orbit 
and away from the pre-SN progenitor orbital velocity. 
For the uniform-direction Vr prior, the kick is directed 
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Table 1. Confidence ranges on progenitor parameters 



Vr prior Median likelihood value 68% confidence interval 95% confidence interval 
Vfc m.2,i <5 Vfe m2,i S Vk m2.i <5 
(km/s) (Mq) (°) (km/s) (Mq) (°) (km/s) (1%) (°) 

J0737-3039A/B 

Gaussian 165 1.80 12.0 80-305 1.50-2.40 3.0-24.5 45-1005 1.37-4.00 1.0-102.5 

Uniform-in-cos 60 1.45 3.5 40-80 1.37-1.55 1.5-5.5 20-140 1.37-1.80 0.5-11.0 

B1534+12 

Gaussian 235 2.45 24.5 200-270 2.00-3.35 20.5-28.5 175-305 1.60-3.90 16.5-32.5 



nearly randomly in the plane perpendicular to the pre-SN 
progenitor orbital velocity. 

The only likely radial velocities are either slightly nega- 
tive or else positive. This reflects the fact that large negative 
velocities (relative to our current reference frame) imply that 
the system must have been born in the outer reaches of the 
Galaxy, but there are far fewer potential birth sites at large 
Galactic radii. 

Low, even very low (< 10°), misalignment angles are 
predicted between A's spin axis and the orbital angular 
momentum, in excellent agreem ent with t he obse rved lack 
of pr ofile variations llManchcste r et all 1200 5: Kra mer et alJ 

The immediate B progenitor was probably less than 
about 2Mq. Since the pre-SN orbital period was compa- 
rable to the current 0.1-day value, the low mass is likely due 
to si gnificant rn a ss los s acc ompanying the orbi tal shrink- 
age. iDewi et alJ J2002l) and iDewi fc Po3 J2003l) trace the 
histories of He-star/NS binaries, showing that suitable pa- 
rameters can be obtained from a range of starting points. 
A difficulty is that the lowest-mass systems are expected 
to undergo common-envelope evolution and spiral-in, re- 
sulting in extremely tightly-b ound or even merged systems 
JPewi fc Polj (jiool), but see llvanova et al.1 J2003M . Dewi 
& Pols note, however, that the relative time-scales of the 
spiral-in and supernova explosion are not well known, and 
an explosion might occur before the spiral-in is complete. 

The low-mass B progenitor is consistent with either an 
electron-capture colla pse of the ON eMg core of a low-mass 
He-star ( Podsiadlowsk i~ t al.ll20o'5ll or the collapse of a low- 
mass iron core. Either might occur sufficiently fast to pre- 
vent the d evelopment of large asymmetries, resulting i n the 
small kick llPfahl et al.l2002t |Podsiadlowski erai]l2004) . Es- 
timates of the time-scales for the onset of the two types of 
explosions are similar (De wi fc Pol3l2003l) and appear too 
long relative to the spiral-in time. Better modeling of the 
two types of explosions may ultimately allow us to distin- 
guish between the two on the basis of the required short 
time- scale. 

PSR B1534-(-12 presents an interesting contrast. Our 
derived pdfs are fully compatible with the parameter space 
allowed in TDS05; in particular the kick velocity is likely 
in the range 200-270 km/s while the progenitor companion 
mass was in the range 2.0-3.35 A/0. We report the confi- 
dence intervals only for the gaussian (dispersion 200 km/s) 
Vr prior in Table the uniform-in-cos prior gives simi- 
lar results. We note that the most important constraint 
in this case is the measured misalignment angle (25 ± 4°; 
IStairs et al.l l2004ll . whose only model-dependence is on 



the reasonable assumption tha t the rotating vector model 
iRadhakrishnan fc Cookelll969^ describes the l inear polar- 
ization position-angle swing of this pulsar (c.f. iDewi et alJ 
I2OO51I . Overall, the parameters of the B1534-I-12 progenitor 
fit quite well with a more 'standard' evolutionary picture 
that involves He-star Roche-lobe overfiow but not a late 
spiral-in phase. It is possible that the longer orbital period 
and/or higher mass of the B 1534-1- 12 progenitor permitted a 
larger core to develop before the SN explosion, leading to a 
larger kick; we note that the second- formed NS in B 1534-1- 12 
is 0.1 M0 more mass i ve tha n J0737— 3039B. Also of inter- 
est is that llhm et"ai] ll200d) find PSR B1913-I-16, with the 
largest known second-formed NS mass, cannot be formed in 
population syntheses with small kick velocities. 

The contrast between these systems urges caution in 
discussing expectations based on low kicks. In the last sev- 
eral years, small (a 50 km/s) kicks have been proposed 
to be nearly universal for nascent neutron stars processed 
in relatively close binary systems, due to the early loss of 
the NS progenitor's envelope and subsequent formation of 
a low-mass core which in turn provides a small kick upon 
collapse (e.g.. Podsiad lo wski et al. 2004n). This hypothesis 
has been used to explain the existen ce of long-period, low- 
eccentricity high-mass X-ray binaries jPfahl et all2002r) . the 
reten tion of large numbers of neutron stars in glob ular clus- 
ters JPfahl et alJl2002L iPodsiadlowski et alJl2004^ ■ the ap- 
parent correlation between the orbital ecce ntricity and the 
recycled pulsar spin pe riod in DNS systems llFaulkner et alJ 
l2005l:lDewi et alJl2005^ and the apparent dearth of isolated 
mildly recycled pulsars ('fa iled DNSs') ejecte d from unbind- 
ing second SN explosions iDewi et alj 2005), as well as to 
predict larger numbers of DNS mergers than are inferred 
from the observed set of objects dPfahl et al.ll2003) . Our de- 
rived progenitor parameters for J0737— 3039A/B certainly 
imply that some fraction of DNS progenitors will be low- 
mass and experience small kicks, and hence that each of 
these explanations is plausible at some level. However, the 
much larger kick needed for B1534-I-12, which is also a short- 
period (lOh) binary, indicates that a range of parameters 
must be considered in all of these arguments. Noting the 
additional rough correlation between the orbital eccentric- 
ity and the mass of the second-formed compact object in 
the six systems for which the masses are well-determined 
(PSRs J1756-2251, J0737-3039, J1 141-6545 , B2303+46, 
B1534-I-12 and B1913+16; see also iFaulkneil i2004)), we 
speculate that the magnitude of the kick may depend quite 
sensitively on the size of the collapsing core. This specula- 
tion is not supported by the most recent 2-D numerical sim- 
ulations of hydrodynamic instabilities during SN explosions 
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iScheck et al. llioOfil) . which suggest that core mass and NS 
velocity are not correlated, but more work will be required 
on 3-D simulations and on other kick models. At the same 
time, careful population synthesis, incorporating initial mass 
functions and preferably an improved understanding of the 
supernova vs. orbital-evolution time-scales, will be required 
to determine whether this correlation or the (likely related) 
spin-period-eccentricity relation can truly be reproduced by 
a mass-dependent kick. 

As pointed out bv IPodsiadlowski et al.l ||2QQ^, the low 
space velocity of J0737— 3039 may provide a further probe 
of its long-term evolutionary history. A competing model to 
the 'standard' DNS evolution ary model out l ined in Section 1 
is the 'double-core' scenario llBrownlll995l : iBethe fc Brownl 
llQQST) in which the He cores of two main-sequence stars of 
nearly equal mass undergo a spiral-in through the envelopes 
of both stars simultaneously. Thus at the time of the first 
SN, the secondary star is a (low-mass) Ife star rather than a 
massive main-sequence star, and the systemic velocity after 
the first supernova explosion would be expected to be higher 
than in the standard model. The observed low space velocity 
therefore makes the double-core scenario less likely. If SN 
kicks are in fact lower in He stars processed in binaries, 
this argument is slightly weakened; however we note that 
under the mass-dependent kick hypothesis, since A's mass is 
close to that of the companion in B1534+12, its kick velocity 
might well have been large. 

Finally, our low-velocity progenitor for J0737— 3039A/B 
may have implications for the nature of short GRBs. In the 
popular DNS inspiral models (IPacvnskilll98(t) . a poor corre- 
lation between active star formation and GRB activity can 
be a natural consequence of the time delay between birth 
and inspiral and the space velocity of the post-supernova 
binary. The few detected counterparts to short GRBs place 
these objects in the outskirts of a broad v ariety of galaxies 
witho ut active current star formation ('e.g:.. |Prochaska et alJ 
l200ti) . Allowing for large DNS kicks, this is consistent with 
the produc tion of the progenitors d uring normal galactic star 
formation jBelczvnski et afll2006ri . However, if smaller DNS 
kicks are assumed, then it may also be necessary to produce 
a significant fraction of the progenitors in globular clusters 
in the galaxy haloes (iGrindlav et alJl2006l) . Our low velocity 
for J0737— 3039A/B tends to support the latter scenario, but 
also suggests that similar systems could be a source of GRBs 
near the Eart h, with possible implications for our terrestrial 
environment jThorsettHligllThomas et alJl2005^ . 
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